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Section  1 


INTRODUCTION 

This  user's  manual  provides  instructions  to  initiate,  run,  and  terminate  the 
External  Burning  Propulsion  Analysis  (EBPA)  computer  program.  This  program 
was  developed  by  the  McDonnell  Douglas  Astronautics  Company  under  soonsorship 
of  the  Air  Force  Rocket  Propulsion  Laboratory,  Edwards  Air  Force  Base, 
California.  The  program  was  developed  for  operation  on  the  CDC  64D0/7600  and 
CYBER  174  computing  systems,  but  its  machine  dependency  has  been  minimized. 
The  EBPA  program  is  self-contained  and  requires  no  communication  with  outside 
programs.  In  addition  to  the  program  tape,  a  single  chemistry  data  tape  is 
required  for  its  operation. 


I 


Section  2 


INPUTS 


The  EBPA  program  is  essentially  composed  of  two  computational  modules:  a 
Fuel  Rich  Jet  Interaction  (FRJI)  subprogram  and  the  Base  Burning  Modeling 
(BBMODL)  subprogram.  The  input  data  for  the  EBPA  program  will  be  described, 
first  for  those  inputs  required  to  execute  the  FRJI  portion  of  the  program, 
which  will  define  the  flow  field  at  the  base  plane,  and  second,  for  the  inputs 
exclusively  required  for  the  injectant  annulus  surrounding  the  base  cavity. 

2 . 1  FRJI  Subprogram  Input 

The  input  data  for  the  FRJI  subprogram  consists  of  vehicle  parameters,  flight 
parameters,  injector  parameters,  chemistry  parameters,  local  conditions,  and 
option  flags.  An  example  of  the  loadsheet  is  shown  in  Figure  1.  The  input 
glossary  is  provided  in  Table  1  which  gives  the  name,  symbol,  location,  and 
definition  of  each  input  variable. 


Figure  2  shows  a  schematic  of  the  biconic  vehicle  configuration.  From  the 

schematic,  the  usage  of  9, ,  9„ ,  X,  .  ,  X  ,  I  ,  and  Z.  .is  straightforward. 
3  12  body  cone  eg  jet 

Since  a  conic  vehicle  is  a  special  case  of  a  biconic  vehicle,  9,  =  9^  for  a 


conical  vehicle.  A  boattailed  vehicle  can  be  approximated  by  setting  9^ 


Flight  parameters  must  be  input  when  the  local  conditions  at  the  injector  are 
to  be  calculated  by  the  computer  program.  When  inputting  free-stream  condi¬ 
tions  M_,  T  ,  Pm,  the  local  condition  flag  LC  in  location  068  must  be  set  to 
1.0  for  the  local  flow  subroutine  to  make  the  injector  computations.  The  pro¬ 
gram  defaults  to  LC  =  1.0.  When  the  free-stream  conditions  (M  .,  Tv,  and  Pj 
are  input,  the  vehicle  angle  of  attack  a  and  roll  anqle  :  must  also  be  input. 
It  is  recommended  here  that  zero  be  used  for  both  values.  The  sign  convention 
on  (  and  :  is  such  that  when  lookingforward  with  a  base  view  of  the  vehicle  a 
wind  direction  from  a  counterclockwise  direction  is  posi ti ve.  Wi til  the  injector 
located  on  top  (9.  =  0),  a  positive  rt  will  put  the  injector  on  the  windward 
ray,  and  a  negative  .  with  9 j  =  0  will  put  the  injector  on  the  leeward  ray. 
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EXTERNAL  BURNING  PROPULSION  ANALYSIS 


LOAOSHEET  1 
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Figure  1  EBPA  Sample  Loadsheet  1 
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If  the  local  conditions  on  the  vehicle  are  known,  it  is  not  necessary  to  input 
Ma>,  T^,  P^,  a,  or  <j>.  Instead,  the  local  conditions  P^,  T^,  and  are  input, 
and  flag  LC  set  equal  to  0.0.  It  is  not  necessary  to  input  as  this  is 
always  calculated  internally.  Note  that  in  either  type  input  (free  stream 
or  local),  the  specific  heat  ratio  for  air  y  must  always  be  input. 

The  type  of  fuel  KINDF  and  type  of  oxidizer  KINDO  in  the  chemistry  parameters 
can  be  selected  from  the  list  of  fuels  in  data  block  FRJIB,  as  shown  in  the 
Appendix.  The  parameters  ETAR  and  CCRIT  are  determined  from  experimental  data. 
The  number  of  zones  Nz  can  vary  from  1  to  10,  but  3  are  generally  sufficient. 
The  oxidizer  to  fuel  ratio  0/F  is  the  fractional  ratio  of  oxidizer  to  fuel.  In 
the  event  that  the  fuel  type  already  contains  the  oxidizer  as  in  some  solid 
propellants,  input  the  same  fuel  type  number  for  both  KINDF  and  KINDO  and  set 
0/F  equal  to  1.0. 

The  option  flags  consist  of  two  print  flags  PFft  and  PFg,  and  the  local  condi¬ 
tion  flag  LC.  All  are  given  in  the  input  glossary  in  Table  1. 

2 . 2  BBMODL  Subprogram  Inputs 

These  inputs  consist  of  geometry  parameters,  annulus  confi'guration  control 
variables,  and  base  bleed  parameters.  An  example  of  the  loadsheet  is  shown  in 
Figure  3.  The  input  glossary  is  included  in  Table  1. 

Geometry  inDuts  include  the  number  of  jet  nozzles  around  the  vehicle  periphery, 
streamline  divergence  angle  at  the  base,  mixing  angle,  and  inner  flow  to  base 
plane  radius  ratio.  Annulus  configurations  are  determined  by  divergence 
streamline  angles  downstream  of  the  base,  cavity  wake  half  angle  and  inner  flow 
half  angle.  These  inputs  are  used  to  indicate  if  system  performance  is  to  be 
determined  or  to  provide  results  from  tabulated  values  of  the  angular  measure¬ 
ments.  The  table  option  can  be  used  to  make  estimates  of  the  region  of 
convergence  of  base  pressures.  The  print  flag  PFp  should  be  1 .  or  2.  if  the 
table  option  is  exercised.  Base  bleed  inputs  include  flow  rate,  base  bleed 
area,  and  either  bleed  flow  pressure  drop  or  bleed  flow  veloci ty ,  one  of  which 
must  be  input  if  base  bleed  is  specified  (N0). 
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Table  1  (Page  1  of  5) 
EBPA  INPUT  GLOSSARY 


Name 

Symbol 

Location 

Definition 

THETA1 

91 

001 

Cone  half-angle.  When  the  vehicle  is  an 
elliptic  cone,  it  is  the  half-angle  along 
the  minor  axis.  For  a  biconic  vehicle  it 
is  the  forebody  half-angle  (deg). 

THETA2 

92 

002 

Body  half-angle.  When  the  vehicle  is  a 
cone  this  angle  should  be  input  equal  to 

THETA1 .  For  a  biconic  vehicle  it 
represents  the  aft  body  half  angle  along 
the  major  axis  (deg). 

XBODY 

^body 

003 

Overall  vehicle  length.  In  the  case  of  a 
cone  or  elliptic  cone  vehicle,  the  overall 
length  is  also  equal  to  XCONE  (ft). 

XCONE 

Y 

cone 

004 

The  vehicle  forebody  length.  The  length 
of  the  conic  section  for  a  biconic  vehicle. 

In  the  case  of  a  cone  or  elliptic  cone 
shaped  vehicle  this  is  also  equal  to  the 
overall  body  length  XBODY  (ft). 

MINF 

M 

03 

005 

Free-stream  Mach  number. 

ALT 

ALT 

006 

Flight  altitude.  This  variable  is  not  used 
in  any  computations.  It  is  for  printout 
only  (ft). 

PINF 

P 

oo 

007 

Free-stream  pressure  (psf). 

TINF 

T 

oo 

008 

Free-stream  temperature  (°R). 

ALFA 

a 

009 

Angle  of  attack  (deg) . 

PHI 

PHI 

010 

Wind  angle  or  roll  angle  (deg). 

ICG 

rcg 

Oil 

Vehicle  center  of  gravity  (ft). 

GAMA 

Y 

012 

Specific  heat  ratio  of  air. 

POJET 

P. 

J 

013 

Thruster  chamber  pressure  (psi). 

AOJET 

Aj 

014 

p 

Thruster  nozzle  throat  area  (in  ). 

ERATIO 

^ratio 

015 

Ratio  of  nozzle  exit  area  to  nozzle  throat 
area. 
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Table  1  (Page  2  of  5) 
EBPA  INPUT  GLOSSARY 


Name 

Symbol 

Location 

Defini tion 

ALFJ 

a . 

J 

016 

Forward  inclination  of  jet,  measured  from 
the  normal  to  the  vehicle  surface  (deg). 

ZJET 

Zjet 

017 

Thruster  site  distance  measured  from  the 
vehicle  base  (in). 

BETA 

R 

018 

Half-angle  of  nozzle  exit  cone. 

KINDF 

KINDF 

019 

Type  of  fuel.  Select  from  one  of  the 
compositions  listed  in  data  block  FRJIB. 

KINDO 

KINDO 

020 

Type  of  oxidizer.  Select  from  one  of  the 
compositions  listed  in  the  data  block 

FRJIB. 

ETAR 

nR 

021 

Indicates  frozen  or  equilibrium  flow. 
O-Frozen  flow 

1-Equilibrium  flow 

ETAR  may  also  be  a  fraction  of  equilibrium. 

PC  TR 

PCTr 

022 

Modifies  the  expression  correlating  the 
upstream  part  of  the  jet  shock.  Optimum 
values  range  between  ,7and  .9.  The  default 
is  .  8 . 

MO 

ml 

024 

Local  Mach  number. 

PPO 

PL 

025 

Local  pressure  (psf). 

PTO 

tl 

026 

Local  temperature  (°R). 

PVO 

VL 

027 

Local  velocity  (ft/sec). 

CCRIT 

Ccrit 

029 

The  percentage  of  fuel  to  be  considered  in 
the  Gaussian  distribution. 

DNI 

DNI 

030 

Input  '  ot  requi red. 

NZONES 

Nz 

032 

Number  of  mixing  zones. 

IPNT 

pfa 

039 

Print  flag 

0-Level  one  printout 

1 - Level  two  printout 

2- Level  three  printout 
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EBPA  INPUT  GLOSSARY 


Name 

Symbol 

Locati on 

Definition 

IPNT 

pfb 

040 

Chemistry  print  flag 

0-No  chemistry  output 

1 -Chemistry  printout  at  each  computa- 
tational  station 

PROG 

*c 

050 

Input  not  required. 

OOVERF 

O/F 

064 

Oxidizer  over  fuel  ratio  of  the  injectant. 

THE TAJ 

9j 

065 

Angular  location  of  the  thruster  relative 
to  the  vehicle  centerline  as  viewed  from 
the  rear  looking  forward.  A  positive 
direction  is  counterclockwise  (deg). 

REFRAD 

RE 

066 

Input  not  required. 

RN 

RN 

067 

Vehicle  nose  radius  (ft).  Defaults  to 

0.001  ft  if  not  input. 

LOCFLG 

LC 

068 

Local  condition  flag 

0-Local  conditions  are  input 

1-Free-stream  conditions  are  input  and 
local  conditions  are  computed  in  sub¬ 
routine  LOCAL 

NNOZZ 

^nozz 

096 

The  number  of  individual  nozzles  around 
the  vehicle  peri phery ,  near  the  base. 

PHIBD 

^B 

097 

The  divergence  angle,  measured  with  res¬ 
pect  to  the  body  surface,  of  streamline 
entering  the  annulus  external  boundary  at 
the  base  pi ane  (deg) . 

SIGBD 

°B 

098 

Mixing  maximum-slope  half  angle  for  air 
and  fuel  injectant  (deg). 

Cl 

CB 

0^9 

Fractional  outer  flow  expansion  angle  from 
base.  If  zero  or  not  input,  no  expansion 
occurs;  full  expansion  if  1.0  is  input. 

PRTFLG 

PFP 

101 

Print  flag  for  base  flow  computations. 

0.,  Summary  Report  only 
=  1.,  Intermediate  iterations,  report. 

2.,  All  iterations,  report. 

RWF 

RWf 

102 

Ratio  of  experimental  reattachment  radius 
to  base  plane  ^adius.  If  not  input,  the 
value  is  computed  by  an  equation  in  the  code. 

Table  1  (Page  4  of  5) 
EBPA  INPUT  GLOSSARY 


Name _ Symbol _ Location 


Definition 


DELTAO,  6 
I  o 


DELTA0o  6 

2 

DELTAF,  6 

1  rl 

DELTAF2  6f 

DELINC,  6 

1  A, 


DELINC0  <5 

2  A. 


103 


104 

105 

106 
107 


108 


Initial  estimate  of  the  divergence  angle, 
measured  with  respect  to  the  body  surface, 
of  streamline  entering  the  annulus  exter¬ 
nal  boundary  at  the  turning  point  for  annu¬ 
lus  1.  The  estimate  is  used  to  start  the 
streamline  divergence  interations.  The  de¬ 
fault  value  is  0.  Optimally,  a  table  of 
equi -distant  intervals  can  be  generated 
(see  DELTNC,  below)  using  DELTAO,  as  the 
first  table  entry  (deg). 


Same  as  6 


for  annulus  2. 


Final  table  entry  for  the  streamline  diver¬ 
gence  angle  for  annulus  1  if  the  table 
option  is  used  (deg). 


Same  as  for  annulus  2. 

T1 

Increment  for  streamline  divergence  angle 
for  annulus  1  if  the  table  option  is  used. 
The  table  option  is  invoked  if  DELINC,  is 
input  non-zero.  The  equi-distant  table  en¬ 
tries  will  consist  of  DELTAO, ,  incremented 
in  steps  of  DELINC, ,  up  to  but  not  ex¬ 
ceeding  DELTAF-j. 

Same  as  5.  for  annulus  2. 

A1 


THETW0,  0,  109 

1  w„ 


THETW0o  0  110 

2  w„ 


Initial  estimate  of  the  cavity  half  angle, 
measured  with  respect  to  the  vehicle  axis 
for  annulus  1.  The  estimate  is  used  to 
start  the  solution  process  of  determining 
system  performance.  A  value  greater  than 
zero  should  be  input.  Optionally,  a  table 
of  equi-distant  intervals  can  be  generated 
(see  THWINC, ,  below)  using  THETW0,  as  the 
first  table  entry  (deg). 

Initial  estimate  of  the  innerflow  half 
angle,  measured  with  respect  to  the  vehicle 
axis,  for  annulus  2;  used  to  start  itera¬ 
tions  for  the  inner  flow  calculations.  Its 
value  must  not  be  zero  if  flow  for  the  2nd 
annulus  is  to  be  generated  (deg). 
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Table  1  (Page  5  of  5) 


EBPA  INPUT  GLOSSARY 


Name 

Symbo 1 

Location 

Definition 

THETWF-J 

*+-” 

CD 

111 

Final  table  entry  for  cavity  half  angle 
for  annulus  1  if  the  table  option  is  used 
(deg). 

thetwf2 

Qw,. 

'2 

112 

Reserved . 

THWINC] 

6w 

A1 

113 

Increment  for  cavity  half  angle  for  annulus 

1  if  the  table  option  is  used.  The  table 
option  is  invoked  if  THWINC]  is  input  non¬ 
zero.  The  equi-distant  table  entries  will 
consist  of  THETWO] >  incremented  in  steps 
of  THWINC] ,  up  to  but  not  exceeding  THETWF], 

The  default  for  THWINC]  is  0. 

thwinc2 

0w. 

A2 

114 

Reserved . 

DELTAR 

6R 

115 

Step  in  annulus  inner  boundary  radius  at 
inner  flow  geometry  transi tion;  al lows  non¬ 
zero  inner  flowrate  in  annulus  2  when  base 
bleed  is  zero  (ft). 

MDOTB 

mB 

116 

Base  bleed  flow  rate  (lb/sec). 

VMATCH 

VM 

117 

Defines  the  point  on  the  inner  shear  layer 
at  which  the  recirculation  velocity  profile 
is  matched  in  velocity  and  gradient.  Input 
presently  not  required. 

AOB 

AOB 

118 

p 

Total  area  of  base  bleed  apertures  (ft  ). 

DELPBI 

apb 

119 

Bleed  flow  pressure  drop.  If  specified 
(>0),  bleed  flow  velocity,  VR,  is  calcu¬ 
lated  (lb/ft2). 

VBI 

VB 

120 

Bleed  flow  velocity.  If  specified  (>0), 
bleed  flow  pressure  drop,  APg,  is  calcu¬ 
lated  (ft/sec). 
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EXTERNAL  BURNING  PROPULSION  ANALYSIS 
LOADSHEET  2 


KEYPUNCH:  INPUT  A-2 


Figure  3  EBPA  Sample  Loadsheet  2 


Section  3 


PROCESSING  FLOW 


3.1  Executive  Routine  EBPA 

The  EBPA  computer  program  uses  an  executive  routine  to  control  two  individual 
subprograms  FRJI  and  BBMODL  as  shown  in  Figure  4.  In  each  subprogram,  control 
is  maintained  by  a  driver  routine  which  establishes  the  proper  sequencing  and 
calling  logic  for  subordinate  subroutines.  An  input  data  editing  subroutine 
(INPUTA)  is  first  called  which  reviews  the  data  and  sets  up  the  data  in  terms 
of  cases,  reference  runs,  and  basic  decks  and  stores  the  edited  data  on  disc. 

A  single  case  is  then  input  to  core,  and  the  two  subprograms  are  sequentially 
executed  based  on  the  input  data.  When  the  computations  are  complete,  control 
is  returned  to  the  INPUT  routine  for  the  next  case  of  data.  If  there  is  none, 
program  execution  is  terminated. 

3.2  Fuel  Rich  Jet  Interaction  Subprogram  FRJI 

A  flow  diagram  for  the  FRJI  subprogram  is  given  in  Figure  5.  The  subprogram 
is  activated  by  the  executive  control  routine.  The  input  data  arelisted,  sub¬ 
program  constants  are  computed,  and  variables  are  initialized.  The  local  flow 
conditions  immediately  upstream  of  the  injection  nozzle  are  calculated  by  sub¬ 
routine  LOCOND  for  given  free-stream  conditions,  angle  of  attack  and  bank 
angle.  The  conservation  equations  for  energy,  mass,  and  momentum  in  the  up¬ 
stream  region  are  solved  in  the  subroutine  UPCON. 

If  the  solutions  of  the  equations  converge,  a  summary  of  the  upstream  region 
is  printed;  otherwise,  subprogram  control  is  returned  to  the  executive  routine 
for  new  input  data. 

With  the  computations  for  the  upstream  region  completed,  the  program  proceeds 
to  compute  the  downstream  region.  The  subroutine  controlling  the  downstream 
calculations  is  DNCON.  A  test  is  first  made  to  determine  if  the  downstream 
interval  being  calculated  is  the  final  one.  The  outer  flow  field  is  then 
calculated  and  also  the  cross-sectional  area  of  the  inner  flow.  At  this  point, 
the  conservation  equations  for  the  downstream  region  are  solved  by  iteration 
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;  i 


EBPA 


y 


INPUTA 


Input  all  cases  of 
data  and  save  on 
Unit  1 


*1 

JL. 


INPUT 


Read  into  core  one 
case  of  data 


.  ’  .  frji 

Perform  Fuel  Rich 
Jet  Interaction 
Computations 

_ l _ 

_ BBMODL _ 

Perform  Base  Burning 
Analysis 


Figure  4  Flow  Diagram  of  Executive  Control  Subroutine  EBPA 
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I 


Figure  5  Flow  Diagram  for  Subprogram  FRJI 
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with  chemistry.  When  convergence  of  the  equations  for  the  inner  flow  is 
achieved,  a  comparison  is  made  of  the  pressures  at  the  combustion  boundary. 

If  the  pressures  across  the  boundary  are  not  equal,  an  adjustment  is  made  in 
the  boundary  slope  and  the  process  of  computing  the  outer  flow  field  and  the 
inner  region  is  repeated.  When  convergence  in  the  outer  loop  is  achieved,  a 
comparison  is  made  of  the  pressures  at  the  combustion  boundary.  If  the  pres¬ 
sures  across  the  boundary  are  not  equal,  an  adjustment  is  made  in  the  boundary 
slope  and  the  process  of  computing  the  outer  flow  field  and  the  inner  region 
is  repeated.  When  convergence  in  the  outer  loop  is  achieved,  control  is 
returned  to  the  subprogram  FRJI,  and  the  subroutine  OUTER  is  entered  at  OUTERS 
and  the  variables  in  the  outer  flow  field  are  dimensional i zed .  Forces  are 
computed  for  the  downstream  segment  and  a  summary  printout  is  made.  If  the 
downstream  segment  is  not  the  last  one  on  the  vehicle,  the  above  process, 
starting  with  solution  of  the  downstream  conservation  equations,  is  repeated. 
After  the  final  segment  has  been  processed,  a  tabulation  of  the  flow  field 
summary  is  produced.  The  subprogram  is  then  returned  to  the  executive  con¬ 
trol  routine. 


Base  Burning  Modeling  Subprogram  BBMODL 


This  subprogram  computes  the  flow  in  the  separated  region  behind  the  vehicle 
base  plane.  The  flow  conditions  generated  at  the  base  by  the  FRJI  subprogram, 
and  user  inputs,  constitute  the  inputs  to  the  subprogram.  A  flow  diagram  of 
BBMODL  is  given  in  Figure  6. 


The  inputs  from  loadsheet  2  are  printed  and,  if  no  case  termination  error  has 
occurred  upstream  of  the  base,  processed  to  determine  the  method  of  flow  gener¬ 
ation.  A  case  termination  error  causes  control  to  return  to  the  executive 
routine.  Inner  and  outer  flow,  flow  rates  and  thermochemistry  at  the  base 
plane  are  then  saved,  and  working  variables  initialized. 


A  value  of  the  cavity  half-angle  for  the  first  annulus  is  calculated.  A  value 
for  the  streamline  divergence  angle  at  the  turn  for  the  first  annulus  is  then 
computed.  The  annulus  number  is  set,  initial  and  constant  flow  parameters, 
constant  geometry,  and  outer  flow  expansion  computations  are  made  for  the 
annulus  and,  if  requested,  printed.  Chemistry  parameters  are  also  retrieved, 
if  required,  for  the  first  annulus. 
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Figure  6  Flow  Diagram  for  Subprogram  BBMODL 
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A  test  is  then  made  to  determine  if  flow  calculations  are  to  be  made  for  annu¬ 
lus  1  or  2.  If  computations  are  for  the  first  annul  us ,  annul  us  geometry ,  momen¬ 
tum,  forces ,  base  pressure,  and  the  flow  field  for  the  current  cavity  half  angle 
9w-|  and  streamline  divergence  angle  6,  are  generated  and,  if  requested, 
printed.  If  no  case  termination  error  occurred  or  a  new  6-j  is  not  required  as 
a  result  of  low  temperature  or  geometry  inconsistency,  the  inner  flow  just 
before  reattachment  is  calculated.  If  a  case  termination  error  occurs,  the 
subprogram  is  returned  to  the  executive  control  routine  for  anew  case  of  data. 
If  a  new  6-j  is  required,  6-|  is  advanced  and  computations  repeated  using  the 
new  6,.  If  the  calculated  base  pressure  Pg-|  (or  P.+.|)  differs  relatively  from 
the  outer  bounding  pressure  Pout^  downstream  of  the  turn  in  flowby  less  than  a 
pre-assigned  tolerance,  Pg-j  is  considered  converged  for  the  current  9w^  and  ^  . 
Calculations  then  proceed  to  the  second  annulus. 

If  a  pressure  solution  does  not  yet  exist,  and  the  maximum  number  of  attempts 
advancing  6-j  has  not  been  used,  a  new  6-j  is  computed  and  the  flow  calculations 
repeated.  If  the  maximum  number  of  attempts  in  6^  has  occurred,  a  diagnostic 
is  printed,  a  new  Qw-j  is  determined  and  all  calculations  for  the  annulus  are 
then  repeated.  If  the  maximum  number  of  computations  for  a  Qw-j  occurs,  a  diag¬ 
nostic  is  printed  and  control  is  returned  to  the  executive  routine  for  a  new 
case  of  i  nput  data  . 


The  second  annulus  is  entered  on  convergence  of  Pg-j  (annulus  1).  Chemistry 
parameters  of  annulus  1  are  saved.  Initial  and  constant  flow  parameters, 
constant  geometry,  and  outer  flow  constants  for  annulus  2  are  aenerated 
and,  if  required,  printed.  An  inner  flow  half  angle  9w^  for  the  annulus  is 
estimated.  A  value  for  the  streamline  divergence  angle  6^  at  the  flow  turn  is 
then  computed.  The  annulus  geometry,  conservation,  flow  field  and  a  new  9w^ 
from  the  geometry  are  computed.  Output  is  printed  if  requested.  If  no  case 
termination  error  occurred,  the  pressure  Pgj  at  the  critical  point  is  tested 
for  convergence.  A  case  termination  error  causes  control  to  return  to  the 
executive  routine  for  a  new  case  of  input  data.  If  there  is  P$2  convergence. 


a  test  is  made  to  determine  the  change  in  the  previous  and  present  9w^  at  the 
P$2  convergence.  If  the  change  is  within  tolerance,  the  relative  system  error 


between  the  inner  flow  pressure  Pw  at  reattachment  and 


is  computed  and 
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tested.  If  the  error  is  within  tolerance,  the  system  solution  has  been 
achieved,  a  Summary  Report  of  the  results  is  output,  and  return  is  then  made 
to  the  executive  routine  for  a  new  case  of  input  data. 

If  Pf32  has  not  converged  for  the  current  0W£  and  ^  a  new  6?  f°und  and 
flow  computations  repeated.  If  the  maximum  number  of  computations  for  has 
occurred,  a  diagnostic  is  printed,  9w^  is  updated  from  the  latest  computation 
and  the  annulus  2  computation  process  is  repeated  with  the  current  Gw^.  If 
there  is  convergence  but  change  in  the  previous  and  present  9w^  at  the 
P^2  convergence  is  not  within  tolerance,  9w£  is  updated  to  the  present  value 
and  annulus  computations  repeated.  When  the  maximum  number  of  0w^  updates 
occurs,  a  diagnostic  is  printed  and  the  system  error  calculated.  If  the 
relative  system  error  in  Pq^  and  P^  is  not  within  tolerance,  control  is 
returned  to  annulus  1,  a  new  9w-|  for  the  first  annulus  is  determined  from 
the  system  error,  and  the  entire  process  of  the  above  calculations  is 
repeated.  The  repetitions  continue  until  the  system  error  is  within  toler¬ 
ance  or  the  maximum  number  of  computations  for  9w-|  has  occurred. 
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Section  4 


INSTRUCTIONS  FOR  USE 


4.1  Program  Activation  and  Control 

The  following  is  an  example  of  the  control  card  sequence  of  instructions  and 
deck  setup  required  for  activation  and  execution  of  the  EBPA  program. 


To  execute  the  program  on  the  CDC  6600  from  tape  with  source  card  images: 
CC1 


J0B,T400,P2,MT2.  (Accounting) 

LABEL, EBPA,R,D=HY,L=EBPAPS,VSN=Reel  #. 
REQUEST, TAPEW,HY,VSN=Ree1  #. 
REWIND(EBPA.TAPEW) 

FTN( I=EBPA,R=3,PL=50000) 

RETURN(EBPA). 

MAP, ON. 

SEGLOAD. 

LDSET,PRESET=NGING. 

LGO. 

7-9-8  Card 

SEGLOADer  directives 
See  below 

7-8-9  Card 
Data  Cards 
7-8-9  Card 
6- 7-8-9  Card 


Source  Card  images  file 
Chemistry  Tape 

Compile  with  full  cross 
reference  and  page  limit. 

Prints  full  loader  map. 
Invoke  program  segmentation 
Set  core  to  negative 
infinites. 

Execute  program. 

Multi -punch  end  of  record 
card. 


Multi-punch  end  of  job  card 


The  default  file  names  for  the  input,  output  and  chemistry  tape  files  are, 
respectively,  INPUT,  OUTPUT,  and  TAPEW.  If  the  program  input  data  resided 
on  file  INPDAT,  the  program  output  written  on  file  OUTDAT,  and  file  CDAT 
contained  the  chemistry  tape,  the  LGO  card  would  become: 

LG0( INPDAT, OUTOAT, CDAT)  . 

The  SEGLOADER  directives  for  the  CDC  7600,  CYBER  174  and  CDC  6t  >  are 
listed  below.  For  the  CYBER  and  6600,  XTOI  =  and  IOCON.  are  not  required. 
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r 


CC1 

ecu 

CC21 

EBPaT 

TREE 

=  BPA«<r«Jl,B8'10DLi  INPUT*) 

EBPA 

SBPA 

INCLUDE  lfeC0D6-,SN0ODb.lFLTQUT,,SNDFrL,(E0F*XTOl. 

INCLUDE  FLTlN«,iHcOh»iINPC»iA«AKER»»DUTB* 

SBPa 

ESPa 

-'StDBAlr- 

GLOBAL 

global 

SlLL-iGVm - 

CUE'UiCHE'42,c*  =  u3»LH£>'4  *C0NST2 

D'iC0N1.0NCU':2,3sUATAlEsncV 

-wi — 
EBPa 
EBpA 

GLOBAL- 

global 

GLOBAL 

TTCTOBtF£TGbt I'FLGWjMUUTSiFUfcL'  - 

GIRLP«H6An«INPgT,UfeT,jPE,N»LnCALS- 

MAlN.’1AlNli*Al''2|MAsl,MAs2»HAS3»MISA,MAs5 

espa 

EBPA 

global 

GLOBAL 

global 

-frrAT.-T'fXa,  .1TX3,  taRK"«3  ,  IOC  PIP  ,'-I  1 1"» 

0 J fESL I OuT£0  2. cu fEP3, 0UTE«4 1  OUTERS iOUTEB6iOUTEP7»:uTERa 
3z£R0(s30OY,s!eATAiSUMSAF'iTEMPTF;<  *RAS*»uPCDn1 

— cspt“ 
ebPa 
ebPa 

global- 

GLOBAL 

GLOBAL 

■Z5«c,zsiGijr 

C3NE,DCOSNe.FRsTRH#THR06E 

a».I0.jl0CON.,i0.BUF.#ALLlh.HLlSTJJJFLAG  ....  ... 

- Ei<nr~ 

- - - 

4.2  Preparation  of  Input  Data 

The  EBPA  program  includes  a  general  purpose  data  input  editor  subroutine  which 
converts  data  from  punched  cards  to  machine-internal  formats  acceptable  to 
the  program.  The  subroutine  provides  a  method  of  loading  numeric  and  alpha¬ 
betic  information  into  core,  of  partitioning  the  data  into  groupings,  and  of 
minimizing  the  high  rate  of  redundancy  in  data  that  has  to  be  read  into  the 
program.  This  Basic  Deck-Reference  Run-Case  hierarchy  of  data  groupings  is 
discussed  below.  The  following  material  discusses  the  data  card  formats  per¬ 
mitted  by  the  program. 

4.2.1  Card  Formats 
Format  1  -  Job  Title  Cards 

This  input  format  is  used  for  loading  alphameric  data.  Figure  7.  Column  1 
contains  an  8-3  (’’=")  punch.  Columns  2-5  contain  the  location  within  the 
array  in  which  the  first  data  word  is  to  be  stored.  Columns  6-53  are  used 
for  loading  alphameric  information,  eight  words  per  card.  These  eight 
words  in  columns  6-11,  12-17,  18-23,  24-29,  30-35,  36-41,  42-47,  and  48-53 
will  be  stored  in  the  relative  locations  LOC  through  L0C+7,  respectively. 
Columns  54-62  should  not  be  used. 
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Figure  7.  Job  Title  Card 

The  standard  80-column  FORTRAN  coding  sheet  may  be  used  in  lieu  of  a  load- 
sheet.  Column  1  must  contain  an  equal  sign  ("=");  columns  2-5  contain  the 
location  within  the  data  array  of  the  first  word  on  the  card;  columns  6-53 
contain  the  eight  six-character  words  to  be  loaded;  and  columns  63-72  carry 
the  BD,  RR,  and  CASE  numbers.  Columns  54-62  should  not  be  used. 

Format  2  -  Numeric  Data 

This  format  is  designed  for  entry  of  real  number  values.  A  card  guide  for 
entering  the  values  is  shown  in  Figure  8.  Column  1  must  contain  a  12-5 
("E")  punch.  Columns  2-61  are  used  for  entering  real  number  data,  five 
fields  per  card.  Columns  2-4,  14-16,  26-28,  38-40,  and  50-52  contain  the 
location  within  the  array  in  which  the  data  is  to  be  stored.  Columns  5-13, 
17-25,  29-37,  41-49,  and  53-61  contain  the  values  of  the  elements. 

In  the  loadsheet  entry  for  this  format,  the  Column  1  space  must  contain  an 
"E".  The  QUAN  column  is  for  the  programmer’s  reference  only,  and  is  not 
to  be  punched.  The  LOC  column  is  used  to  identify  the  location  of  the  data 
within  the  one-dimensional  data  array.  The  values  of  the  elements  are 
entered  in  the  VALUE  column,  up  to  nine  characters,  including  the  sign, 
decimal  point,  and  E  type  exponent,  if  used. 
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Figure  8  Numeric  Data  Cara 

The  permissible  variations  in  the  way  these  data  may  be  entered  are: 

A.  Plus  signs  may  either  be  punched  or  omitted,  but  minus  signs  must 
be  punched. 

B.  For  numbers  using  E  conversion,  (i.e.,  2.5E01),  the  start  of  the 
exponent  field  must  be  indicated  by  an  "E",  or  by  a  plus  or  minus 
sign--the  start  of  the  field  cannot  be  indicated  by  a  blank. 

C.  The  decimal  point  must  be  punched  to  insure  that  the  value  will  be 
interpreted  correctly. 

D.  The  exponent,  if  any,  must  be  right  justified  in  the  field. 

E.  The  INTEGER  value  in  the  LOC  field  must  be  right  adjusted. 

The  following  applies  to  aJJ_  input  data  cards: 

Cards  with  a  zero  or  a  blank  in  column  1,  and  cards  with  any  punch 
other  than  a  zero  or  a  blank  in  column  62  will  be  ignored.  Columns 
63-65  contain  the  Basic  Deck  number;  columns  66-67  contain  the 
Reference  Run  number;  and  columns  68-70  contain  the  CASE  number. 
Columns  71-76  should  remain  blank.  Columns  77-80  may  be  used  for 
the  job  sequence  number  or  other  identification.  All  data  cards 
must  be  in  sort,  major  to  minor,  columns  63-70. 
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4.2.2  Basic  Deck-Reference  Run-Case  Data  Input  Capabilities 

A  proper  definition  of  terminology  is  essential  to  a  discussion  of  data  input 
editor  capabilities.  Many  misconceptions  on  this  subject  are  che  result  of 
using  terminology  which  is  loosely  defined.  To  avoid  such  confusion,  the 
terminology  for  this  discussion  is  defined  below.  The  user  is  advised  to 
refer  to  this  list  of  basic  definitions  from  time  to  time. 

Case  Study  -  A  computerized  simulation  of  some  process  for  a  defined 
set  of  parameters. 

Run  •  The  groups  of  data  required  for  one  or  more  case  studies,  and 
which  are  collectively  submitted  at  one  time  to  be  processed  by  a 
computer  program. 

Set  of  Data  -  The  data  that  are  required  for  a  case  study. 

Data  Element  -  One  of  the  items  of  data  in  a  set  of  data. 

BASIC  DECK  -  REFERENCE  RUN  -  CASE  Hierarchy  -  The  hierarchical 
organization  of  data  elements  used  in  a  case  study.  This  hierarchy  is 
based  on  the  degree  of  constancy  each  particular  data  element  exhibits 
from  one  case  study  to  the  next. 

Data  Level  of  a  data  element  -  The  place  in  the  BASIC  DECK  -  REFER¬ 
ENCE  RUN  -  CASE  hierarchy  at  which  that  data  element  is  defined  for 
a  particular  case  study.  If  a  data  element  is  defined  in  a  BASIC  DECK 
(BD),  a  REFERENCE  RUN  (RR),  or  a  CASE,  its  data  level  is  respectively 
the  BD-level,  the  RR-level,  or  the  CASE-level. 

BD  Data  -  The  collec.  . ■  of  data  elements  which  is  defined  at  the  BD- 
level,  and  which  contributes  to  the  typical  set  of  data  being  discussed 
(versus  all  data  elements  in  the  run  which  are  defined  at  the  BD-level). 
The  term  is  sometimes  shortened  to  "BD."  BD  data  have  the  highest  degree 
of  constancy  from  one  case  study  to  the  next. 

RR  Data  -  The  collection  of  data  elements  which  is  defined  at  the  RR- 
level,  and  which  contributes  to  the  typical  set  of  data  being  discussed 
(versus  all  data  elements  in  the  run  which  are  defined  at  the  RR-level). 
The  term  is  sometimes  shortened  to  "RR."  RR  data  have  some  lesser 
(than  BD)  degree  of  constancy  from  one  case  study  to  the  next. 


29 


CASE  Data  -  The  collection  of  data  elements  which  is  defined  at  the 
CASE-level,  and  which  contributes  to  the  typical  set  of  data  being 
discussed  (versus  all  data  elements  in  the  run  which  are  defined  at 

the  CASE-level).  The  term  is  sometimes  shortened  to  "CASE"  or  "CS." 

CASE  data  are  unique  to  a  case  study. 

Data  Origin  -  The  location  in  computer  memory  where  the  first  element 

of  a  group  of  related  data  elements  is  stored. 

Part  Case  Data  -  The  collection  of  data  elements  which  results  when  a 
set  of  data  is  partitioned  into  several  smaller  groups.  This  term  is 
sometimes  shortened  to  "PC."  PC  data  can  occur  at  any  level  of  the 
BD-RR-CASE  hierarchy. 

Domain  of  an  RR  -  Consists  of  the  grouping  of  CASE  data  for  all  case 
studies  which  use  that  RR  to  define  data  elements. 

Domain  of  a  BD  -  Consists  of  the  groupings  of  RR  data  and  CASE  data  for 
all  case  studies  which  use  the  BD  to  define  data  elements. 

Basic  Deck,  Reference  Run,  and  Case 

BD,  RR,  and  CASE  capability  is  designed  to  facilitate  handling  of  data  input 
from  punched  cards  where  there  is  a  high  degree  of  redundancy. 

Concept  of  Reference  Run  and  Case 

Because  of  the  nature  of  computer  applications  often  encountered,  experience 
has  shown  that  a  typical  program's  user  will  read  a  quantity  of  data  which 
describes  a  basic  model;  the  user  will  then  change  the  values  of  a  few 
parameters  and  generate  "case  studies"  for  different  variations  from  that 
basic  model.  The  purpose  of  the  RR  and  CASE  concept  is  to  provide  such  a 
user  the  facility  to  refer  back  to  a  particular  data  set  used  earlier  (the 
basic  model)  and  to  instruct  the  computer  to  use  that  data  set  aqain,  with 
certain  specified  changes.  The  typical  data  set  will  thus  contain  data 
elements  from  two  sources;  RR  data  consist  of  those  data  elements  of  a  set  of 
data  which  stay  the  same  (i.e.,  which  are  referenced)  from  one  case  study 
of  a  basic  model  to  the  next;  CASE  data  consist  of  those  data  elements  of  a 
set  of  data  which  are  unique  to  a  particular  case  study. 
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Example  of  Reference  Run  and  Case 

Suppose  one  wanted  to  simulate  the  flight  of  a  missile  over  a  gi ven  trajectory , 
varying  several  parameters  which  define  the  missile's  shape  and  calculating 
the  resultant  effect  on  the  missile's  surface  temperatures.  In  this  instance, 
the  trajectory  data  would  ideally  be  defined  as  the  RR  data;  the  missile  shape 
parameters  would  constitute  CASE  data.  One  set  of  data  representing  one 
simulated  flight  cf  the  missile  over  the  trajectory  would  be  defined  by  the 
RR  data  and  a  particular  CASE;  a  second  set  of  data  representing  a  simulated 
flight  of  a  modified  missile  over  the  same  trajectory  would  be  defined  by 
the  RR  data,  but  with  a  different  CASE.  The  collection  of  the  RR  data  and 
the  individual  CASEs  (which  comprise  that  RR's  domain)  might  constitute  the 
run. 

Use  of  Reference  Run  and  Case  Data 

When  using  the  INPUTA  subroutine,  the  user  specifies  which  data  items  are 
RR  data  and  which  are  CASE  data.  By  careful  selection  as  to  what  he 
includes  in  the  RR  category,  the  user  can  minimize  the  chore  of  data  trans- 
scription  onto  loadsheets.  The  user  must  follow  certain  conventions  when 
using  RR  and  CASE  capability.  Failure  to  follow  these  conventions  or  failure 
to  have  a  clear  understanding  of  the  process  may  result  in  the  computer 
program  either  indicating  a  data  error  or  interpreting  the  data  in  an  undesi red 
manner. 

A.  Both  RR  and  CASE  data  are  identified  by  numbers.  These  numbers 
are  assigned  by  the  user  and  are  physically  punched  in  the  cards 
which  contain  the  corresponding  data  elements. 

1.  The  user  defines  a  set  of  RR  data  by  assigning  a  two-digit 
nonzero  RR  number  (e.g.,  "01,"  "27,"  "83").  The  CASE  number 
field  for  that  RR  data  must  be  defined  as  "000."  These  two 
numbers  will  be  punched  in  all  data  cards  which  define  that  RR. 

2.  The  user  defines  CASE  data  by  assigning  a  three-digit,  nonzero 
CASE  number  (e.g.,  "001,"  "270,"  "508").  He  must  also 
specify  the  applicable  RR  number.  These  two  numbers  will  be 
punched  in  all  data  cards  containing  the  data  elements  which 
define  that  CASE. 
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B.  If  the  user  wants  to  study  more  than  one  basic  model  in  a  run,  he 
can  define  as  many  as  99  sets  of  RR  data,  such  that  each  set  has  a 
unique  RR  number.  CASE  number  assignments  need  be  unique  only 
within  the  domain  of  a  given  RR.  Thus,  one  run  could  contain  the 
two  unrelated  sets  of  data  identified  as  "RR  1  CASE  5"  and  "RR 
CASE  5). 

C.  An  RR  number  of  "00"  exists,  but  cannot  be  used  to  define  RR  data. 
Such  a  designation  is  punched  on  CASE  data  cards  when  a  run  (or 
part  of  a  run)  is  being  made  which  does  not  require  RR  capability, 
or  in  instances  when  BD  is  being  used  (to  be  discussed  later). 

D.  Both  the  RR  and  CASE  fields  on  the  punched  data  cards  must  be 
completely  punched  with  numerics  (the  digits  0  through  9).  When 
computer  operations  keypunches  the  data,  blanks  are  replaced  by 
zeros . 

E.  The  load  sheet  forms  define  the  RR  number  (occupying  two  card 
columns  for  its  two  digits)  and  the  CASE  number  (occupying 

three  card  columns  for  its  three  digits)  in  two  adjacent  fields;  the 
RR  number  field  is  to  the  left  of,  and  hence  higher-ordered  than, 
the  CASE  number  field. 

F.  The  data  cards  for  a  run  must  be  ordered  so  that  successive  cards 
present  RR  and  CASE  numbers  which,  if  read  together  as  one 
number,  define  a  monotonical ly  increasing  sequence.  RR  data  must 
be  defined  in  the  sequence  of  data  cards  prior  to  any  CASE  data  in 
its  domain.  All  CASE  data  cards  that  refer  to  a  given  RR  number 
must  be  contiguous  and  ordered,  and  must  adjoin  those  data  cards 
which  define  the  RR. 

G.  The  computer  program  will  not  accept  RR  data  by  itself.  Program 
execution  using  a  set  of  data  is  initiated  by  the  reading  of  CASE 
data.  The  minimum  CASE  consist  of  one  punched  card  with  at  least 
one  data  item  punched  on  that  card  in  the  left-most  data  field. 

H.  When  the  user  is  assigning  his  data  to  the  RR  and  CASE  categories, 
he  should  remember  that  a  set  of  data  is  effectively  generated  by 
reading  the  RR  data  into  memory,  then  overlaying  it  with  those  data 
items  defined  in  the  CASE  data.  By  this  overlay  process,  data 
items  not  defined  in  the  RR  data  are  assigned  a  value  if  defined 

in  the  CASE  data;  data  items  defined  in  the  RR  data  are  assigned 
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the  substitute  value  as  (and  if)  redefined  in  the  CASE  data;  data 
items  not  defined  in  either  the  RR  or  CASE  data  are  assigned  a 
value  of  zero.  Successive  sets  of  data  in  that  RR's  domain  are 
generated  in  the  same  manner,  the  CASE  data  always  overlaying  the 
original  RR  data. 

I.  As  a  corollary  to  item  G,  it  is  seen  that  in  a  given  RR's  domain, 
one  CASE  is  independent  of  another.  A  given  CASE  in  that  domain 
can  temporarily  (i.e.,  for  that  set  of  data)  change  items  defined 
in  the  RR. 

Concept  of  Basic  Deck 

The  purpose  of  BD  capability  is  to  provide  the  RR-CASE  user  with  an  extension 
of  that  facility.  BD  data  has  a  relationship  to  RR  data  that  is  similar  to 
the  relationship  of  RR  data  to  CASE  data.  The  BD-RR-CASE  concept  is  useful 
for  programs  whose  data  structure  is  complicated  in  the  degree  of  its 
relative  constancy  from  one  case  study  to  the  next.  In  such  instances,  a 
typical  set  of  data  will  contain  data  elements  from  three  sources:  BD  data, 
RR  data,  and  CASE  data.  As  a  general  rule,  BD  data  consist  of  those  data 
elements  of  a  set  of  data  which  have  the  highest  degree  of  constancy  from 
one  case  study  to  the  next.  RR  data  consist  of  those  data  elements  of  a 
set  of  data  which  have  some  lesser  degree  of  constancy  from  one  case  study 
to  the  next. 

Example  of  Basic  Deck 

Suppose  one  wanted  to  simulate  the  flight  of  three  different  missiles  over  a 
given  trajectory,  varying  several  parameters  defining  each  missile's  shape, 
and  calculating  the  resultant  effect  on  each  missile's  surface  temperature. 

In  this  instance,  the  trajectory  data  would  ideally  be  defined  as  the  BD  data; 
each  missile's  basic  shape  parameters  would  constitute  an  RR.  Variations 
in  these  shape  parameters  would  constitute  CASE  data.  One  set  of  data 
representing  one  simulated  flight  of  the  first  basic  missile  shape  over  the 
trajectory  would  be  defined  by  the  BD  data,  the  RR  data  for  that  missile 
shape,  and  the  CASE  data  defining  one  variant  of  the  shape  parameters.  A 
second  set  of  data  representing  a  simulated  flight  of  that  same  basic  missile, 
but  with  slight  modifications,  over  that  same  trajectory,  would  be  defined  by 
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the  BD  data,  the  same  RR  data,  but  with  a  different  CASE.  A  third  set  of 

data  representing  a  simulated  flight  of  the  second  basic  missile  shape  over 

the  same  trajectory  would  be  defined  by  the  same  BD  data,  but  with  a  second 

set  of  RR  data,  and  one  of  the  CASEs  in  that  second  RR's  domain.  The 

collection  of  BD,  RR's  and  CASE'S  might  constitute  the  run. 

Use  of  Basic  Deck  Data 

When  INPUTA  is  used  to  provide  a  BD  capability,  the  following  rules  apply. 

The  rules  defined  for  RR  and  CASE  remain  applicable,  unless  specifically 
noted  otherwise.  The  user  must  follow  certain  conventions  when  using  BD 
capability.  These  conventions  are  simply  extensions  of  those  defined  in  the 
RR  and  CASE  discussion,  with  similar  consequences  if  not  practiced. 

A.  BD  data  are  identified  by  number.  BD  numbers  are  assigned  by  the 
user  and  are  physically  punched  into  the  cards  which  contain  the 
corresponding  data  elements. 

1.  The  user  defines  a  set  of  BD  data  by  assigning  a  three-digit 
BD  number  (e.g.,  "021,"  "146,"  "938").  The  RR  and  CASE 
number  fields  for  that  BD  data  must  be  respectively  defined 
"00"  and  "000."  These  three  numbers  will  be  punched  in  all 
data  cards  which  define  the  BD. 

2.  In  defining  RR  data,  the  user  must  specify  the  applicable  BD 
number.  That  BD  number  will  be  punched  in  all  data  cards 
which  define  that  RR.  Similarly,  in  defining  CASE  data,  the 
applicable  BD  and  RR  numbers  must  be  punched  in  all  data 
cards  which  define  that  CASE. 

B.  As  many  as  999  unique  sets  of  BD  data  can  be  defined  in  one  run. 

RR  number  assignments  need  be  unique  only  within  the  domain  of 

a  BD.  Thus,  one  run  could  contain  the  two  unreleased  sets  of  data 
identified  as  "BD  1  RR  5  CASE  6"  and  "BD  2  RR  5  CASE  6." 

C.  A  BD  number  of  "000"  exists  but  cannot  be  used  to  define  BD  data. 
Such  a  designation  is  punched  on  RR  and  CASE  data  cards  when  a  run 
(or  part  of  a  run)  is  being  made  which  does  not  require  BD  capabi 1 i ty. 
Similarly,  an  RR  number  of  "00"  can  be  used  when  a  CASE  need  only 
refer  to  a  BD. 
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D.  The  BD  field  on  punched  data  cards  must  be  completely  punched 
with  numerics.  When  computer  operations  keypunches  the  data, 
blanks  are  replaced  by  zeros. 

E.  The  loadsheet  forms  which  allow  the  BD  capability  define  the  BD 
number  field  (occupying  three  columns  for  its  three  digits)  to 
the  left  of  and  adjacent  to  the  RR  number  field. 

F.  The  data  cards  for  a  run  must  be  ordered  so  that  successive  cards 
present  BD,  RR,  and  CASE  numbers  which,  if  read  together  as  one 
number,  define  a  monotonical ly  increasing  sequence.  BD  data 
must  be  defined  in  the  sequence  of  data  cards  prior  to  any  RR  or 
CASE  data  in  the  BD  domain.  RR  and  CASE  data  cards  that  refer 

to  a  given  BD  number  must  be  contiguous  and  ordered,  and  must 
adjoin  those  data  cards  which  define  the  BD. 

G.  The  computer  program  wil 1  not  accept  BD  or  BD/RR  data  by  itself. 
Program  execution  using  a  set  of  data  is  initiated  by  the  reading 
of  CASE  data.  The  minimum  CASE  consists  of  one  punched  card 
with  at  least  one  data  item  punched  on  that  card  in  the  left-most 
data  field. 

H.  When  the  user  is  assigning  his  data  to  the  BD,  RR  and  CASE 
categories,  he  should  remember  that  a  set  of  data  is  effectively 
generated  by:  (1)  reading  the  BD  data  into  memory,  (2)  overlaying 
it  with  those  data  items  defined  in  the  RR  data  (if  any),  and  (3) 
overlaying  this  "intermediate  result"  with  those  data  items  defined 
in  the  CASE  data.  (The  overlaying  process  is  the  logical  extension 
of  that  described  in  RR  and  CASE.  RR  has  the  same  relationship  to 
BD  that  CASE  has  to  RR.)  Successive  sets  of  data  in  that  BD's 
domain  are  generated  in  the  same  manner,  the  RR  and  CASE  data 
always  successively  overlaying  the  original  BD  data. 

I.  As  a  corollary  to  item  G,  it  is  seen  that  RR's  in  a  given  BD's  domain 
are  independent  of  each  other.  A  given  RR  in  a  BD's  domain  can 
temporarily  (i.e.,  for  all  CASE'S  in  its  domain)  change  items 
defined  in  the  BD.  A  given  CASE  in  an  RR's  domain  can  temporarily 
(i.e.,  for  that  set  of  data)  change  items  defined  in  its  associated 
RR  or  BD. 
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Examples  of  BD-RR-CASE 

Some  examples  of  BD-RR-CASE  number  sequencing  of  data  cards  are  as  follows: 


BD 

RR 

CASE 

000 

00 

001 

1 

set  of  case  data  only 

000 

00 

001 

000 

00 

004 

2 

sets  of  case  data  only 

000 

00 

006 

000 

01 

000 

000 

01 

002 

1 

RR  with  2  data  cases 

000 

01 

003 

001 

00 

000 

001 

00 

003 

1 

BD  with  2  data  cases 

001 

00 

005 

002 

00 

000 

002 

04 

000 

1 

BD  with  1  RR  and  1  Case 

002 

04 

005 

005 

00 

000 

005 

03 

000 

005 

03 

008 

005 

03 

009 

1 

BD  with  2  RR's  and  2  cases  for  each  RR 

005 

04 

000 

005 

04 

010 

005 

04 

on 

4.2.3  Chemistry  Tape  Format 

The  chemistry  tape  contains,  primarily,  coefficients  for  computations  of 
enthalpy,  entropy,  and  specific  heat  for  each  of  the  atomic  species  utilized 
by  the  rograin.  The  dato  are  read  in  subroutine  SEARCH  and  the  coefficients 
are  us*  in  subroutine  GIBBS.  The  tape  is  written  in  binary  (unformatted) 
mode  and  tape  density  (bits  per  inch)  is  set  at  generation. 
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Three  records  of  information  are  provided  for  each  species.  The  record 
formats  and  contents  are  described  as  follows. 


Variabl e( s) 

Definition/Usage 

Record  1 

Word  1 

LEVEL 

INTEGER  indicating  coefficient 
data  to  read  if  >  0;  end  of  data 
if  <  0. 

Word  2 

JAN 

INTEGER  specifying  record  length 
for  coefficient  data: 

=1,8  words 
=2,  16  words 

Word  3 

KHASE 

=  0,  flags  specific  counter 

Word  4 

SPECIE: 

Alphanumeric  atomic  symbol  for 
its  species. 

Record  2 

Words  1 

-14 

IEi,  J 

1=1,...,  7 
j=l  ,2 

Atomic  description  INTEGER  array 
providing  for  7  atomic  elements 
3=1 :  number  of  atoms 
j=2:  atomic  number 

Record  3 

Words  1 
or 

-8 

RA,RB ,RC,RD 

Arrays  of  floating  point  co¬ 
efficient  data.  Record  length 

Words  1 

-16 

RA ,  RB ,  RC  ,RD , 
RE,RF,TL  ,TU 

<  10  if  JAN=1  or  16  if  JAN=2. 

The  record  formats  are  repeated  for  each  species. 
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Section  5 


OUTPUT 


Three  levels  of  printed  output  have  been  provided  for  the  EBPA  program.  The 
level  selection  of  printout  is  specified  by  the  PFA  and  PFR  input  flags  in 
locations  039  and  040  of  loadsheet  1  (Figure  1)  andPFpat  location  101  of 
loadsheet  2  (Figure  3).  The  input  values  for  these  flags  are  discussed 
in  Section  2.  Flags  PF^  and  PF^  pertain  to  output  generated  upstream  of  the 
base  plane  by  the  FRJI  portion  of  the  program, and  PFp  to  output  generated 
in  the  annulus  regions  by  the  BBMODL  Subprogram. 

5. 1  Summary  Level  Output 

Tables  2  through  6  provide  the  normal  or  basic  mode  of  output  for  the  pro¬ 
gram.  This  level  of  output  is  obtained  by  a  zero  or  blank  for  the  values  PF^ 
in  location  039  of  loadsheet  1  and  PFp  in  location  101  of  loadsheet  2.  Table  2 
lists  the  user  inputs  from  loadsheet  1  and  consists  of  vehicle,  flight  and 
injector  parameters,  propellant  and  chemical  flow  parameters,  the  determined 
local  and  jet  conditions,  jet  thermochemistry,  and  the  propellant's  chemical 
composition.  Table  3  is  a  summary  of  upstream  region  results  generated  at 
the  completion  of  the  upstream  solution.  Station  data  tabulations  are 
generated  after  all  computations  for  each  station  have  been  made.  A  flow 
field  summary,  Table  4,  tabulates  station  distance,  inner  shock  layer 
radius,  shock  cross-sectional  area,  and  chemical/physical  properties  for 
each  station.  The  tabulations  in  this  table  include  the  last  upstream 
station  and  all  downstream  stations.  The  location  of  the  interface  between 
the  upstream  and  downstream  regions  is  defined  by  the  penetration  height  of 
the  jet. 

Table  5  consists  of  the  inputs  from  loadsheet  2.  These  inputs  are  used 
with  inputs  from  loadsheet  1  and  results  of  the  baseplanefor  calculations 
in  the  base  control  region.  Table  6  is  abbreviated  output  of  the  final 
iteration  in  9  i  for  the  system  solution.  Included  are  the  half  angles, 
divergence  angles,  pressures,  and  the  relative  error  in  the  system  solution. 
Table  7  represents  a  Summary  Report  of  the  EBPA  program. 
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5.2  Intermediate  Level  Output 
5.2.1  Station  Output 

This  level  output  is  produced  upstream  of  the  base  and  is  obtained  by 
setting  flag  PF^  equal  to  1 .  It  includes  all  output  contained  in  Level  1 
plus  additional  tables  of  expanded  information  of  station  data.  The  tables 
appear  for  each  station  in  the  downstream  region  to  the  base  plane. 

The  downstream  station  summary  is  shown  in  Table  8.  It  includes  station 
location,  current  summary  flow  field  data,  and  a  station  flow  field  for  the 
computational  ray  measured  from  the  body  normal  to  the  shock.  This  flow 
field  contains  for  each  point  along  the  ray,  its  Z-distance  and  radius 
measured  in  penetration  heights,  velocity  normalized  to  freestream  velocity, 
pressure,  density  normalized  to  freestream  density,  local  velocity,  tempera¬ 
ture,  stream  function  (PSI)  value,  and  a  measure  of  constant  enthalpy  at  the 
PSI  value.  The  first  flow  field  line  of  output  is  on  the  shock,  and  the 
last  is  at  the  body.  The  force  computations  indicated  in  the  upper  left 
corner  are  deleted. 


5.2.2  Annulus  Output 

This  level  of  output  is  produced  at  and  downstream  of  the  base  plane  and  is 
obtained  by  setting  PFp  equal  to  1 .  Table  9  lists  flow  parameters  and  geometry 
at  the  base  when  computations  are  entered  for  annulus  1  and  conditions  of  the 
annulus  2  entrance  when  computations  are  entered  for  the  2nd  annulus. 

For  annulus  1,  AI,  PI,  MI,  VI,  MDOTI,  MOMI,  and  FI  denote  the 
area,  pressure,  Mach  number,  velocity,  mass  flow  rate,  momentum  and  force, 
respectively,  of  the  annulus  flow  from  each  nozzle  at  the  base  plane.  The 
PIP1  and  MDTIP1  are  pressure  and  flow  rate  at  reattachment  and  are  initialized 
to  the  values  at  the  base.  The  MOLWT,  ROCHEM,  TCHEM,  and  GAMMA  denote 
molecular  weight,  density,  temperature  and  specific  heat  from  thermochemistry. 
The  outer  flow  parameters  are  indicated  by  M ^ »  T1 ,  RH01 ,  VI,  M1P,  PIP, 
RH01P,  and  VIP,  with  the  P-  appended  variables  representing  a  Prandtl -Meyer 
expansion  of  flow  around  the  base.  The  MDOTNI,  MODTNA  and  MDOTNJ  are  the  total 
mass,  air  and  jet  flow  rates  from  the  nozzles  around  the  base  cavity.  The 
mass  fractions  of  fuel  and  oxidizer  are  XMF  c.nd  XMO.  The  AREAB,  RB  and  RI 
denote  the  base  area,  base  radius  and  annulus  thickness  at  the  base  plane. 


AIN  measures  the  cross  sectional  area  between  the  base  and  reattachment 
circumferences.  The  angles  THETAB  and  TUETAW  are  the  body  angles  measured 
with  respect  to  the  vehicle  surface  and  the  cavity  half  angle.  THETAC  is  the 
input  angle  THETA2  discussed  in  the  section  "EBPA  Inputs".  MVRINL  is 
to  be  ignored. 

At  the  entrance  to  annulus  2,  AI,  PI,  MI,  VI,  MDOTI,  MOMI  and  FI  are  updated 
to  contain  the  values  at  reattachment.  MOLWT,  ROCHEM,  TCHEM  and  GAMMA  are 
initially  the  last  results  computed  in  annulus  1.  Ml ,  PI ,  T1 ,  RHD1 ,  VI,  Ml P , 
PIP,  RH01 P ,  and  VIP  become  the  outer  flow  parameters  for  annulus  2.  MDOTNI 
is  also  updated.  MDOTNA,  MDOTNJ ,  XMF  and  XMO  remain  constant  throughout 
iterations.  AREAB,  RB,  RI  and  THETAB  become  values  at  the  reattachment 
cross  sectional  entrance  to  annulus  2.  AIN  is  the  initial  value  of  the  cross 
sectional  area  between  the  circumferences  at  reattachment  and  the  critical 
point.  The  THETAW  is  the  inner  flow  half-angle  and  THETAC  contains  the  current 
value  of  the  cavity  half-angle  of  annulus  1. 

All  of  the  above  tabulations  are  printed,  if  requested,  for  each  new  compu¬ 
tation  of  the  streamline  divergence  angle  of  annulus  1  and  at  the  entrance 
of  computations  for  annulus  2. 

Tables  10  and  11  are  samples  of  compact  iteration  output  for  annulus  1  and 
annulus  2  respectively. 

The  output  for  annulus  1  is  printed  at  convergence  of  the  base  pressure  from 
iterations  of  the  streamline  divergence  angle  DELTA  and  at  the  current 
iterated  cavity  half  angle  THETAW.  For  this  option,  STEP  is  the  number  of 
iterations  in  DELTA  required  for  convergence.  R3  is  the  annulus  radius  at 
reattachment,  AL  and  HL  distances  from  the  base  center  and  along  the  center- 
line  to  reattachment.  LTURN  is  the  distance  to  and  RTURN  the  annulus  radius 
at  the  turning  point.  TURANG  is  the  turning  angle  of  flow.  ALPHA  is  the 
angle  between  an  equivalent  single  ray  and  streamline  upstream.  DAREA  is 
the  area  for  mass  flow  entrained  in  the  outer  sheet  layer.  P I  PI  is  the 
converged  base  pressure  at  reattachment.  POUTU,  POUTD  and  PIN  are  pressures  on 
the  annulus  outer  --upstream  and  downstream  of  the  turning  point--  and  inner 
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boundaries.  PIPM  and  PMI  are  the  pressure  and  Mach  number  obtained  after 
the  flow  is  turned  at  the  base  plane.  EPSP  is  the  relative  error  between  PIP1 
and  POUTD  and  satisfies  the  criterion  for  convergence.  AIP1  is  the  cross 
sectional  area  of  the  flow  at  reattachment.  FIN,  FOUT,  and  FTOTAL  are  inner, 
outer,  and  total  forces;  MDTIP1  is  mass  flowrate  at  reattachment,  DMTOTA  is 
the  amount  of  added  air,  TOTMOM  is  conservation  momentum,  and  DMOMI  is  the 
change  in  mass  flow  up  and  downstream  of  the  turn  in  flow.  DMDTUP  and  DMDTDN 
are  the  mass  flow  rates.  MI P 1 ,  TCHEM,  V I  PI ,  VMAX,  RCHEM  and  ROCHEM  are 
reattachment  Mach  number,  temperature,  velocities,  perfect  gas  law  constant, 
and  density  of  the  flow.  M2,  P2,  V 2,  and  RH02  are  the  outer  flow  parameters 
at  reattachment.  RHOVUP,  RHOVDN,  VOUTUP,  and  VOUTDN  are  the  pV  products 
and  velocities  up-  and  downstream  of  the  turning  point.  PDREL  gives  a 
relative  error  for  successive  iterates  in  PI  PI-  THETWG,  EPSG,  RSTAR,  HLREL, 
which  have  no  meaning  for  annulus  1,  MDTENT  and  MREAT,  not  computed  at  this 
point,  are  all  set  to  zero.  The  variables  PINR,  PAROU,  MUROU,  PAROD, 

MV  REX,  PAR  IN  and  RMOMCK  are  irrelevant. 

In  annulus  2,  the  output  is  printed  at  each  convergence  of  pressure  at  the 
critical  point,  and  until  there  is  convergence  in  the  inner  flow  half  angle. 
The  variables  are  similarly  defined  as  in  annulus  1.  Those  values  defined 
relative  to  the  base  plane  in  annulus  1  are  now  defined  relative  to  the  cross 
sectional  entrance  to  annulus  2;  and  the  values  defined  at  reattachment  are 
now  defined  at  the  critical  point.  THETAW  and  THETWG  are  the  inner  flow 
half  angle  iterates.  The  EPSG  is  the  relative  error  satisfying  angle 
convergence.  MDTENT  and  MREAT  are  entrained  mass  flow  rate  and  Mach  number 
of  the  inner  flow  at  reattachment.  RSTAR  is  the  annulus  thickness  at  the 
critical  point.  HLREL  gives  a  relative  error  for  successive  iterates  in  HL. 

Table  12  lists  computations  for  the  1st  annulus  inner  flow  at  reattachment. 

If  PFp  =  1,  the  output  is  printed  on  convergence  of  the  divergence  angle 
DELTA  at  the  current  cavity  half-angle  THETAW.  If P^p  =  2,  output  is  printed 
at  each  DELTA  iteration.  The  first  line  of  values  are  reprinted  from  Table  10. 
UOXR  is  centerline  flow  velocity  of  the  recirculation  region,  positive 
downstream.  RCAV  is  annulus  inner  flow  radius,  XMDTSL  and  MDTDSL  are  the 
incremental  and  total  momentum  flux  across  the  dividing  streamline. 
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These  values  are  tabulated  at  distances  HLX  along  the  centerline  from  the 
base  plane  to  reattachment.  PIP1N  is  the  pressure  at  reattachment,  MDOTO 
base  bleed  momentum,  MDOTHL  momentum  flux  out  of  the  inner  flow  at  the 
downstream  end,  MDTENT  entrained  mass  flow  rate,  MRE  and  MREAT  the 
reattachment  Mach  numbers  averaged  over  the  exit  flow  of  annulus  1 
and  the  entrance  flow  to  annulus  2,  respectively. 

5.3  Detail  Level  Output 

5.3.1  Diagnostic/Trace  Output 

This  level  of  output  is  called  by  setting  PF^  equal  to  2.  This  level  of 
output  includes  all  the  printout  of  the  other  two  levels  plus  the  inter¬ 
mediate  computations  of  variables  and  arrays  during  the  interaction  processes 
for  the  conservation  equations  in  the  upstream,  downstream,  and  annulus 
regions,  the  thermochemical  computations,  mixing  calculations,  and 
calculations  for  the  beginning  of  the  downstream  region.  This  level  of 
output  is  primarily  for  debugging  and  maintenance  operations  and  examples  are 
not  included  in  this  manual  because  of  the  large  quantity  of  output. 

5.3.2  Chemistry  Data  Printout 

Additional  chemistry  data  may  be  obtained  for  each  station  by  setting  PFD 

D 

equal  to  1.  The  output  consists  of  the  chemical  composition  of  products  for 
each  combustion  zone  at  each  upstream,  downstream,  and  annulus  station. 

5.3.3  Extended  Annulus  Output 

This  output  is  an  extension  of  the  iteration  output  discussed  in  the  section 
"Annulus  Output"  and  is  obtained  by  setting  PFp  equal  to  2.  A  sample  of  the 
output  is  shown  in  Tables  13  and  14.  In  annulus  1,  the  STEP  value  in 
this  option  identifies  a  computation  step.  STEP  0  is  a  first  computation 
for  a  given  THETAW  and  DELTA.  The  1  and  2  STEP'S  are  generated  to  update 
thermochemistry  to  be  consistent  with  the  PIP1.  In  annulus  2,  the  above 
looping  is  repeated  to  achieve  a  consistent  HL  as  reflected  in  the  value 
HLREL.  All  other  variables  have  the  meanings  described  in  the  section 
"Annulus  Output." 

5 . 4  Sample  Case  Loadsheet 

Figures  9,  10  and  1 1  represent  loadsheet  data  entries  to  produce  the  sample  case 
run. 
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Figure  9  Sample  Case  Loadsheet  BD  Entries 


EXTERNAL  BURNING  PROPULSION  ANALYSIS 


LOADSHEET  2 


KEYPUNCH:  INPUT  A-2 


Figure  10  Sample  Case  Loadsheet  RR  Entries 


EXTERNAL  BURNING  PROPULSION  ANALYSIS 
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Figure  11  Sample  Case  Loadsheet  CASE  Entries 


OPERATING  MESSAGES 


This  section  lists  messages  that  are  for  the  user's  information.  The  mess¬ 
ages  indicate  that  although  the  program  terminates  in  normal  fashion,  the 
expected  results  have  been  influenced  by  corrective  actions  within  the  pro¬ 
gram.  The  messages  are  listed  under  the  name  of  the  subroutine  in  which 
they  occur. 

Subroutine  AMTURN 

AMTURN  NON-CONVERGENT  IN  20  ITERATIONS 

Local  conditions  computed  for  biconic  vehicle  may  be  in  error, 
or  Mach  expansion  angle  in  the  annulus  region  is  excessive.  Program  con¬ 
tinues  execution. 

Subroutine  BBMODL 

1.  SOLUTION  NOT  FOUND  FOR  THETAW  n 

2.  SOLUTION  NOT  FOUND  FOR  DELTA  n. 

The  maximum  number  of  iterations  for  the  half  angles  or  the  divergence 
angles  in  the  annulus  region  has  occurred.  The  n  is  the  annulus  number. 
The  program  continues  execution  until  the  message  occurs  for  the  THETAW  1. 
At  that  point,  execution  proceeds  to  the  control  routine  for  a  new 
case  of  input  data. 

Subroutine  GEMTRY 

1.  ***HL  =  XXX  LTURN  =  XXX  AT  DELTA  =  XXX  AND  THETAW  =  XXX*** 

TRY  NEXT  DELTA 

The  distance  HL  to  the  reattachment  (or  critical  point)  is  less  than  the 
distance  LTURN  from  the  base  (or  reattachment)  to  the  turning  point  for 
annulus  1  (or  2).  Program  will  adjust  values  and  continue  execution. 

2.  INITIAL  THETAW  TOO  LARGE  -  TRY  NEW  GUESS 

The  cavity  half  angle  of  annulus  1  will  be  reduced  and  execution  continued. 

3.  ***TRY  NEW  THETAW  THETAW  DELTA  R3  RTURN  HL  L  TURN  RSTAR  AIP1 
XXX  XXX  XXX  XXX  XXX  XXX  XXX  XXX 

The  current  cavity  half -angle  of  either  annulus  has  caused  the  annulus 
radius  and/or  cross  section  area  at  reattachment  or  critical  points  to 


become  negative.  The  XXX's  represent  the  values  listed.  The  program 
will  continue  execution. 

SUBROUTINE  INCRMT 

1.  ***  DELTA  RESET  FROM  XXXX 

2.  ***  THETAW  RESET  FROM  XXXX 

The  Newton-Raphson  method  has  produced  an  increment  of  undesirable 
magnitude.  The  value  is  reset  from  that  indicated  and  execution 
continues.  The  message(s)  are  printed  only  if  t'  e  PFp  flan  of  Loadsheet 
2  is  set  equal  to  2. 

Subroutine  OUTER 

1.  MAXIMUM  ITERATIONS  IN  OUTER 
AT  DOWNSTREAM  STATION  XXX 
SUMMARY  PRINTOUTS  FOLLOW 

Inner  and  outer  solutions  upstream  of  the  base  have  not  converged. 
Summary  printouts  and  totals  will  be  printed  up  to  the  indicated  station 
and  execution  will  proceed  to  the  next  data  case. 

2.  PRESSURE  PERTURBATIONS  BECOME  INSIGNIFICANT 
AT  DOWNSTREAM  STATION  XXX 

SUMMARY  PRINTOUTS  FOLLOW 

Expansion  process  is  such  that  flow  Mach  angle  has  been  achieved  and 
no  further  expansion  can  occur.  Summaries  and  totals  are  printed  and 
the  program  execution  proceeds  to  the  next  data  case. 

Subroutine  PBSOLV 

***  MAX  LOOPS  FOR  HLREL  XXXX  -  WILL  TRY  NEW  DATA 

This  message  applies  to  the  annulus  2  region  and  indicates  inner  loop 
iterations  have  not  produced  a  satisfactory  computation  for  the  reat¬ 
tachment  to  critical  point  distance.  Execution  continues  with  the 
existing  value  of  HL. 

Subroutine  PMOMNT 

***TCHEM  =  XXXX  -  TRY  NEXT  DELTA 
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Temperature  computations  i 
permitted  by  the  program, 
execution  continues.  The 


n  the  annulus  region  have  fallen  below  that 
A  new  divergence  angle  is  determined  and 
message  is  printed  only  when  PFp  =  2. 
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Section  7 
TROUBLESHOOTING 

This  section  provides  information  necessary  for  the  user  to  troubleshoot  the 
subprograms  to  the  degree  needed  to  isolate  and  correct  user  errors.  All 
diagnostic  messages  and  their  explanations  are  listed  under  the  name  of  the 
subroutine  in  which  they  are  printed. 

Subroutine  DEFIQJ 
CANT  FIND  BASIS 

FATAL  ERROR  IN  H099,  EXECUTION  TERMINATED 

An  improper  chemistry  tape  has  been  mounted.  Check  for  correct 
reel  number.  Program  terminates  execution. 

Subroutine  DNCON 

1.  NEWTON  RAPHSON  LOOP  FAILED,  WILL  TRY  AGAIN  VC=XXX.XXXX 

Program  is  attempting  to  estimate  velocity  in  downstream  region.  The 
message  is  information  only  and  the  program  continues. 

2.  EST  VIP1  CONVG  FAIL,  VG,  VC,  VCDR,  VIP1  ,  VMAX,  CONV 
(Printed  values  for  each  variable) 

Program  is  attempting  to  estimate  for  velocity  in  downstream  region. 

A  new  guess  is  made  and  the  program  continues. 

3.  ROCHEM  CONVG  FAILURE,  ROCHEM,  ROBB,  VIP1,  VIP1C 
(Printed  values  for  each  variable) 

Convergence  failure  in  computation  of  velocity  from  continuity.  A  new 
guess  is  made  and  the  program  continues. 

4.  OVER  40  ITERATIONS  IN  DOWNSTREAM 

Iteration  with  chemistry  to  solve  conservation  equations  in  density, 
velocity  and  pressure  have  exceeded  40.  Program  proceeds  to  next  data 


5.  CHOKED  FLOW 

No  solution  exists  at  this  point  of  downstream  iteration.  Program 
assumes  sonic  flow  and  continues. 


Subroutine  MIX 

MIX  WILL  NOT  CONVERGE 

Output  from  routine  unobtainable  in  30  iterations.  Program  requires 
modification  for  rectification.  Discard  data.  Program  will  proceed  to 
next  case  of  data. 

Subroutine  PERFRM 

FATAL  ERROR  IN  PERFRM 

Message  follows  indication  of  error  in  fuel  ingredient  table.  Program 
proceeds  to  next  case. 

Subroutine  SEARCH 

HARK,  NO  COMBUSTION  SPECIES  FOR  (NAME)  REVISE  PEPAUX 

Combustion  species  indicated  could  not  be  found  on  the  chemistry  input 

tape.  Program  proceeds  to  next  case. 

Subroutine  SIBODY 

1.  THIS  NOTE  IS  FROM  THE  RADIUS  SUBPROGRAM,  AND  INDICATES  THE  EXISTENCE 
OF  AN  ERROR  CONDITION 

Applies  to  the  following  messages: 

2.  THIS  SUBROUTINE  CONCLUDES  THAT  THE  SOLUTION  RAIR  IS  LESS  ( GREATER)  THAN  THE 
MINIMUM  (MAXIMUM)  PERMISSIBLE  VALUE  OF  RMIN  ( RMAX) .  THE  MINIMUM 
(MAXIMUM)  PERMISSIBLE  VALUE  WILL  BE  RETURNED  AS  THE  SOLUTION  VALUE. 

Information  message  in  the  computation  of  ip  values  at  the  body.  RAIR 

is  radius  of  curvature  parameter. 

Subroutine  STREAF 

1.  BODY  STREAMLINE  ABOVE  JET  SHOCK 

2.  STREAMLINE  INTERSECTS  SHOCK  BEYOND  Z=1 

For  either  message  there  is  no  solution  for  outer  conditions  with  given 
inner  conditions.  Program  proceeds  to  next  data  case. 

Subroutine  UPCON 

1.  VI  CONVERGENCE  INCOMPLETE 

Information  only.  The  Newton-Raphson  solution  for  upstream  velocity  has 
not  been  achieved.  A  new  guess  is  made  and  another  solution  is  attempted. 
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2.  CHOKE  FLOW 


No  solution  exists  at  this  point  of  upstream  iterations.  Program  assumes 
sonic  flow  and  continues. 

Subroutine  INPUTA 

The  following  error  messages  will  be  printed  if  any  of  these  conditions 
occur  when  input  data  cards  are  loaded. 

1.  CARD  IGNORED 

Printout  of  card;  execution  will  continue 

Either  card  column  1  contains  a  zero  or  blank,  or  card  column  62 
contains  any  punch  other  than  a  zero  or  blank. 

If  any  one  of  the  following  errors  occur,  these  error  messages  will 
be  printed,  and  execution  will  be  terminated. 

2.  PROGRAM  NUMBER  ON  DATA  CARD  INCONSISTENT  WITH  ACTUAL  PROGRAM 
NUMBER  XXXX 

Printout  of  card 

The  program  number  has  been  punched  in  the  data  cards,  and  it 
does  not  agree  exactly  with  the  program  number  referenced  in  the 
CALL  INPUTA  statement.  Leave  columns  73  -  76  blank. 

3.  INVALID  CARD  COLUMN  1  PUNCH 
Printout  of  card 

Column  1  contains  a  punch  other  than  zero,  blank,  5,  6,  7,  D,  E,  I, 

N,  or  =.  EBPA  data  cards  require  the  E  or  =. 

4.  INCORRECT  CASE  IDENTIFICATION  FIRST  CARD 
Printout  of  card 

First  card  contains  either  a  non-zero  RR  and  a  non-zero  CASE,  or  a 
zero  RR  and  a  zero  CASE. 

5.  CARDS  NOT  IN  SEQUENCE. 

Printout  of  card 

RR  or  CASE  numbers  on  data  cards  are  not  in  monotonically  increasing 
order. 

6.  REFERENCE  RUN  NUMBER  CHANGE  BEFORE  NONZERO  CASE. 

Printout  of  card 

RR  number  has  changed  before  a  CASE  data  card  has  been  read. 


8.  NO  CASE  DATA 
Printout  of  card. 

All  data  cards  were  read,  without  encountering  a  non-zero  CASE  number. 

9.  LOCATION  LESS  THAN  OR  EQUAL  TO  ZERO 
Printout  of  card 

LOC  entry  on  a  data  card  was  <0. 

10.  LOCATION  GREATER  THAN  MAXD  -  XXXXX. 

Printout  of  card 

LOC  entry  on  data  card  is  greater  than  maximum  size  specified  in 
the  CALL  INPUTA  statement. 

11.  MAXD  LESS  THAN  4. 

Maximum  size  specified  in  the  CALL  INPUTA  statement  is  <4. 

12.  BASIC  DECK  NUMBER  CHANGE  BEFORE  NON-ZERO  CASE. 

Printout  of  card. 

BD  number  has  changed  before  a  CASE  data  card  has  been  read. 

13.  NO  BASIC  DECK  DATA 
Printout  or  card 

BD  number  had  changed,  and  a  non-zero  CASE  number  appears  on  the 
same  card. 

If  any  of  the  following  messages  occur,  the  user  should  check  card  column 

1  for  an  "E"  on  numeric  data  cards  and  for  an  "="  on  title  heading  cards: 

PART  CASE  PREVIOUSLY  DEFINED  FOR  MATRIX  INPUT 

PART  CASE  PREVIOUSLY  DEFINED  FOR  NON-MATRIX  INPUT 

DIMENSIONS  OF  MATRIX  HAVE  NOT  BEEN  DEFINED 

NUMBER  OF  ROWS  OR  COLUMNS  LESS  THAN  OR  EQUAL  TO  ZERO 

ROW  OR  COLUMN  INDEX  IS  ZERO 

MAXIMUM  DIMENSIONS  OF  MATRIX  (XXXX.XXXX)  EXCEEDED. 

Subroutine  INPUT 

1.  BD  =  XXX,  RR  =  XX,  CASE  =  XXX,  PC  =  XX 

DATA  ARRAY  SIZE  TOO  SMALL  -  XXXXX,  NECESSARY  SIZE  =  XXXXX 
The  storage  space  reserved  by  the  CALL  INPUT  statement  is  not  large 
enough  to  accommodate  the  referenced  array.  The  maximum  input  array  size 
allowed  by  EBPA  is  600.  No  LOC  entry  should  exceed  this  value. 
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2.  NO  ZERO  PART  CASE  DATA 

Columns  71-72  contain  a  non-zero  value.  Part  case  capability  is  not 
provided  by  the  program. 


1 


Subroutine  EQUIL 
ITER  STOP 

Indicates  all  20  iterations  in  equilibrium  computations  were  used. 
Program  continues. 

Subroutine  GEMTRY 

FAILED  TO  RECOMPUTE  THETAW  -  CASE  TERMINATED*** 

The  maximum  number  of  attempts  to  reduce  the  current  cavity  half  angle 
to  achieve  positive  annulus  radius  and/or  cross  sectional  area  has  been 
reached.  Control  is  returned  to  the  executive  routine  for  a  new  case 
of  data.  This  problem  occurs  usually  from  too  large  an  input  value 
for  the  base  bleed  rate. 

Subroutine  GEOM 
DZ  TOO  LARGE 

Projections  of  computational  element  along  the  plate  are  changed  from 
the  ZY  to  the  ZX  plane.  Execution  continues. 

Subroutine  HBAL 

HBAL  FAILURE  I,  HLP,  Wl(4),  HUP,  ETU  FTL 
(Printed  values  for  each  variable) 

Chemistry  could  not  converge  on  current  temperature  for  given  enthalpy. 
Program  continues  with  generated  value.  The  value  may  be  out  of  range 
for  a  particular  zone  but  not  for  the  whole  mix.  If  so,  nonconvergence 
could  cause  case  failure. 

Subroutine  JETPEN 

*PENETRATI0N  HEIGHT  LARGER  THAN  3/4  OF  INJECTION  SITE  DISTANCE. 
HITE=XXX. XXXXXX  FEET* 

Shock  disturbance  exceeds  distance  from  jet  to  base.  Program  proceeds 
to  next  case.  Lower  penetration  height  or  enlarge  body  length  and 
resubmi t  case. 
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Section  8 


ENVIRONMENT 

This  section  identifies  the  hardware  and  support  software  required  for 
execution  of  the  External  Burning  Propulsion  Analysis  Program. 

PROGRAM  LANGUAGE 

The  External  Burning  program  is  written  primarily  in  the  FORTRAN  extended 
programming  language.  Less  than  2  percent  of  the  program  includes  the 
assembler  programming  language  COMPASS. 

FCRTRAN  Compiler 

The  FORTRAN  Extended  (FTN)  compiler  must  be  used  for  the  EBPA  program.  This 
compiler  has  the  capability  to  optimize  program  execution  speed  and  to  reduce 
the  size  of  object  module. 

COMPASS  Assembler 

The  COMPASS  Assembler  is  the  symbolic  programming  language  for  the  Control 
Data  6000,  7000  and  CYBER  170  Series  control  and  peripheral  processors. 

EXISTENT  UTILITIES 

Object  modules  generated  by  the  FORTRAN  and  COMPASS  compilers  are  loaded, 
edited  and  updated  by  existent  utilities. 

LOADER 

The  loader  provides  editing  and  loading  of  object  modules  generated  by  the 
compilers.  In  addition  to  normal  load  functions,  the  loader  provides 
facilities  to  create  overlay  and  segmentation  load  modules. 

OPERATING  ENVIRONMENTS 

The  EBPA  program  has  been  enabled,  with  minimal  modification  of  segmentation 
directives,  to  operate  under  the  following  systems: 


COMPUTER 

CYBER  170/CDC 
CYBER  174 
CDC  6600 


COMPILER 

7600  FTN  4.1+69 
FTN  4.7+485 
FTN  4.7+476 


ASSEMBLER 


OPERATING 

SYSTEM 


COMPASS  3.0-069  NOS  1.4 

COMPASS  3.6-485  NOS  1.3 

COMPASS  3.5-476  NOS/BE  1.3 
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Appendix 

PROPELLANT  FORMULATIONS  CONTAINED  IN  PROGRAM  DATA  BLOCK 

Data  on  a  variety  of  propellants  and  fuels  are  stored  in  data  block  FRJIB  and 
can  be  called  by  user  input  for  the  EBPA  program.  The  propellants  can  be 
selected  from  a  list  which  includes  bipropellants,  momopropellants  and  solid 
propellants.  Air  and  hydrogen  may  also  be  selected  in  addition  to  mixtures  of 
hydrogen  and  nitrogen  at  two  different  temperature  levels.  The  desired 
propellants  can  be  called  by  number  from  the  following  list: 


cal/g 


cal/g 


cal/g 


cal/g 


cal/g 


cal/g 


cal/g 
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8. 


9. 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


AGJ  Solid  Propellant 

C2.  712  H4.  030  °3.  200  Pb0.  003 
Heat  of  formation  -632.  1  cal/g 

Aluminized  Composite  Solid  Propellant 

H0.  0566  C0.  03196  N0.  00249  °0.  00929  A/0.  01112Ci0.  00213 
Heat  of  formation  -196.0  cal/g 

Boronized  Composite  Solid  Propellant 

H0.  07804  B0.  00924  C0.  04617  N0.  00265  °0.  009636  CV  002128 
Heat  of  formation  -216.  3  cal/g 


ARC  UpSTAGE  JIC  Solid  Propellant 

H0.  4799  S.  1873  N0.  05843  °0.  2208  C/0.  05348 


Heat  of  formation  -424.  5 

H2  at  1500  F 

H2 

Heat  of  formation  2767.  6 

H2/N2  at  1500  F  X(H2)  =  0.  30 

H3.0  N7.  0 

Heat  of  formation  286.  0 

H2/N2  at  1500  F  X(H2)  =  0.  667 


H2.0N1.0 

Heat  of  formation 


531.  0 


H2/N2  at  1000  F 


X(H2)  =*0.667 


H2.0N1.0 


cal/g 


cal/g 


cal/g 


cal/g 


Heat  of  formation  338.7  cal/g 

UpSTAGE  EB  Gas  Generator  Solid  Propellant 

C2.  599  H3.  867  °3.  391  N0.  758 

Heat  of  formation  -741.  1  cal/g 


17. 


18. 


Triethylaluminum 

H6.  0  C15.  0  A,1.0 
Heat  of  formation 
Aluminum  borohydride 

A/  B3.  0  H12.  0 
Heat  of  formation 


-291.7  cal/g 


-300.0  cal/g 
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19.  ARC  Solid  Propellant  Arcadene  168 

C. 7969  H3. 1 561  N . 4576  C 1 . 451 1  A11.2973 

Heat  of  formation  -362.62  cal/g 
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